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Abstract 

The nature of Mean Field Solutions to the Equations of Motion of the Chern- 
Simons Landau-Ginsberg (CSLG) description of the Fractional Quantum Hall 
Effect (FQHE) is studied. Beginning with the conventional description of this 
model at some chemical potential fj,Q and magnetic field B corresponding to 
a "special" filling fraction v = lixpjeB = 1/n (n = 1,3,5 • ••) we show that 
a deviation of \i in a finite range around \iq does not change the Mean Field 
solution and thus the mean density of particles in the model. This result 
holds not only for the lowest energy Mean Field solution but for the vortex 
excitations as well. The vortex configurations do not depend on fj, in a finite 
range about (jlq in this model. However when \i — fj,Q < fj,~ r (or \x — fio > ) 
the lowest energy Mean Field solution describes a condensate of vortices (or 
antivortices) . We give numerical examples of vortex and antivortex configu- 
rations and discuss the range of \x and v over which the system of vortices is 
dilute. 



There have been many interesting approaches to understanding the Fractional Quan- 
tum Hall Effect (FQHE) in which a two-dimensional electron gas is subjected to a large 
transverse magnetic field. One such approach which is often called the "Chern-Simons- 
Landau-Ginzberg" (CSLG) approach to the FQHE |L|- |J makes use of the fact that in 
two-dimensions a fermion (such as the electron) can be treated as a boson to which is 
attached an odd number of magnetic flux quanta. This magnetic flux can be the usual 
magnetic flux but more conventionally one introduces an additional Gauge field known as 
the "statistical" gauge field. The electron is then viewed as a boson with an integer number 
of flux quanta of this "statistical magnetic field". This "binding" of the statistical flux to 
the boson is achieved, technically, by the Chern-Simons interaction which is briefly reviewed 
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below. The key observation in this approach to the FQHE is that the most evident plateaus 
occur when the (conventional) magnetic flux per particle is an odd integer: 

2vrp 1 

z/= — = - n= 1,3,5- •• (1) 

eB n 

Thus if we describe each electron as a boson to which is attached an odd number of statistical 
flux quanta, then the mean statistical magnetic field will cancel the applied magnetic field 
precisely when Eq. [I] is satisfied. The resulting theory is that of an interacting system of 
bosons with no net applied field which will Bose condense. The resulting phenomenology 
has been studied in detail in Refs. fj]- ||. 

There has also been considerable work on describing the Fractional Quantum Hall system 
when the density (or the magnetic field) is perturbed away from the special filling fraction 
described by Eq. |1]. In the CSLG description vortices occur in the system to accommodate 
these fluctuations in the density or in the magnetic field. It has been shown that such 
vortices do exist and they can be found by perturbing the ground state of the system at the 
special filling fraction by adding (or removing) particles to the system. In this way single 
vortex solutions can be found jl], []5J. These vortices are in fact the Laughlin quasi-particles 
and quasi-holes . When they occur in pairs they are neutral excitations analogous to the 
rotons of 3-D superfluids. 

Vortices also form when the external magnetic field is changed with the number of par- 
ticles kept fixed. In this case the perturbation is a constant, uniform shift of the field rather 
than a local density perturbation. Thus, in the limit of infinite volume, a lattice of vortices 
covering the entire space is expected to form just as would happen if the mean particle den- 
sity would be perturbed from that of Eq. (|I|). The conventional way to analyze this system 
is to begin at the special filling fraction (i.e. to consider first the case in which the density 
has been modified so that 2itp/eB = 1/n) and to accommodate the excess density by means 
of these vortex excitations. In this paper we shall take a slightly different approach. The 
electron density is described by means of a chemical potential p. When the magnetic field 
is changed away from the special filling fraction we shall try to find the Mean Field ground 
state of the resulting system by looking for a solution to the resulting equations of motion. 
Since, in this case, the mean Statistical Field does not cancel the applied field there will be 
no uniform mean field solution and the Mean Field equations are extremely "frustrated". 
We shall discuss the properties of solutions to this system of equations analytically and then 
present some numerical solutions. 

Our main goal in this paper is to display some novel features of the CSLG description 
of the FQH system. To this end we adopt the simplified version of the model which is 
described in Ref. Q which allows us to show these features more clearly. We begin with a 
gas of bosons represented by a (nonrelativistic) scalar field i/j with a mass m. The fact that 
these bosons are at some finite density p will be implemented by considering the system at 
some chemical potential p. The system will also have a Statistical Gauge Field whose 
main purpose is to attach an integer number of statistical magnetic flux quanta to each 
boson thus allowing them to describe an electron gas. Our main simplification relative to 
the actual physical situation is to imagine that these bosons have a hard core repulsive 
interaction (oc \ip\ A ) instead of a Coulomb interaction. These bosons are also subjected to a 
fixed external magnetic field B which is described by an electromagnetic potential A^. The 
Lagrangian for this system is given (in units for which h — c — 1) by: 



2 



£ = -i^(do - t(eA - ga )ip + — |(V - z(eA - #a))Vf 

! ' 2 . £fluX a^a x + ^ij) 2 -^ (2) 
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We shall see in a moment that in order for if) to describe electrons n must be an odd integer. 
The third term in this Lagrangian is the Chern-Simons term which is responsible for the 
fact that Statistical Magnetic Flux is attached to each particle. The equation of motion for 
the field a , 

J- b = M (3) 

implies that to each particle (with / d 2 xif)^if) = 1) is attached a Statistical magnetic flux 
g& = 2im. Thus if n is an odd integer each boson has an odd number of Statistical flux 



quanta attached to it so that it has Fermi statistics [H 
In the special case when the filling fraction 

v = ^ = ~ (4) 

eB n K J 

(where p = if>^if> is the density of particles) the average effect of the Statistical field precisely 
cancels that of the externally applied magnetic field i.e. g < b >= eB so that a gauge can 
be chosen in which, on the average, g < Oj >= eAi. This, of course, corresponds to the 
special filling fraction and leads to the CSLG description of the Fractional Quantum Hall 
Effect. 

We begin our discussion by writing down the equations of motion which result from 
varying if), a and in the Lagrangian @. Since is the externally applied Gauge 
Potential due to a spatially constant magnetic field we choose Aq = 0. 

D 2 

i(d o +iga o )ip + 1 ^ip-\(ip^)ip + tMip = (5) 



(^D^-(D^) = ^e ij (-doa j + d j ao) = 7 f-e ij eP (7) 
2m 2-kh 2nn 

where Di = d{ — ieAi + igai is the covariant derivative and ej is the statistical electric field. 
The first equation (|5|) is the generalized (nonlinear) Shroedinger equation for this system. 
The second equation (||) implements the constraint that there are n flux quanta per particle 
as discussed above. The third equation which also results from the Chern-Simons term in 
the action relates the statistical electric field to the particle current. 

The above equations of motion are operator equations for the Quantum Fields if> and 
dj. The technique which we use to study these Quantum Equations (see Ref. Jl|]- 0) is to 
first find a Mean Field Solution to the above equations. The properties of the Quantum 
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system can then be studied by expanding about this mean field. We thus begin by solving 
the equations of motion (^-0) as if they were classical equations. When 

£ = — (8) 

this set of equations has the constant solution \ip\ 2 = j, ej = 0, and gb = eB. If either 
or B is changed so that Eq. (|8|) is not satisfied it would seem difficult to find solutions 
which simultaneously minimize the potential and produce no net field. We now analyze the 
equations in this case in a more systematic manner. 

We begin with several important observations about these equations which are general- 
izations of theorems proven in Ref |TI| in a somewhat different context. The first observation 
is that unless fi/X — eB/(2iTn) (which corresponds to the special filling fraction) there are no 
nonzero covariantly constant solutions (solutions with D^ip = 0) to this system of equations. 
The proof is as follows: first write ip = ^e tQ . Now D^ip = implies that = and that 
dpQ = eA^ — ga^ so that the combination eA — ga is a pure gauge. It follows that eB = gb. 
Now Eqs. (p|) and (|6]) with D^ip = imply that either ip = or that = eb/(2im) so 
that /i/A = eB/(2nn). This completes the proof that nonzero covariantly constant solutions 
only exist when n/\ = eB/(2im). 

In addition to the lack of solutions with D^ip = there are also no nonzero covariantly 
static solutions (i.e. solutions with D ip = 0) unless fi/X = eB/{2im). The proof again 
begins by writing ip = £e* n . D ip = then implies that <9o£ = and that d^Q = eA Q — ga^ = 
—ga . Using this and the third equation of motion Eq. (|7|) as well as the fact that OqA^ = 
we find that 

d {djQ + gaj - eA,) oc e jk {d k Q + ga k - eA k ) £ 2 (9) 

Since £ is independent of time we can identify three possibilities. The first option is £ = in 
which case tp = 0. A second possibility is that djQ + gaj — eAj = in which case gb = eB. 
In this case Eq. (|6|) (the Chern-Simons constraint) implies that \ip\ 2 = eB/2nn which is 
the condition for the special filling fraction. This further implies that £ is spatially constant 
which, from Eq. (|5p forces \ip\ 2 = /i/A. The third possibility for satisfying Eq. ([|), which 
is necessary when we have a nonzero xp with /i/A ^ eB/(2im) is that djfl + gaj — eAj is a 
nonzero oscillating (sinusoidal) function of time. It follows that gb — eB and thus gb is an 
oscillating function of time which then implies that £ oscillates with time which contradicts 
the equation dot; = 0. It follows that except at special filling fractions there is no nonzero 
solution which is covariantly static. 

Despite the fact that no covariantly static solution exists it is possible to find nonzero 
time independent solutions (with d^ip = 0) but with a 7^ 0. As is shown (in a somewhat 
different context) in Ref. ( [O) this can be done by extremizing the energy functional 



H = J d 2 xH = J d 2 x 



(10) 



with Di = di — ieAi + igcii but subject to to the constraint Eq. (Q) that 



Vxa = yH 2 (11) 
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In other words when Eq. (|T0|) is varied with respect to ip and a* subject to the constraint 
( PI) one obtains precisely the equations of motion with d ip = but with a ^ 0. 

The term proportional to a$ appears when eij is varied with respect to ip as required by the 
constraint Eq. ([□]). a will be a solution to Eq. (|7|) with <9 cti = (which can always be 
arranged since ip is time independent). 

Our job is now to find a configuration of ip and which minimizes the energy (|10D 
subject to the constraint ( JTT] ) for various values of the parameters p and B. The simplest 
case which has already been discussed is the case in which p/X = eB/(2im). In this case the 
configuration with p — p,/X both minimizes the potential and forces eb = gB which admits 
a configuration with Ditp = 0. In this situation the density p satisfies v = 27rp/eB = 1/n 
which is the special filling fraction. The case of interest to us however occurs when either 
B or p is modified so that 

In this case there is a problem. If we force the potential to be at its minimum \ip\ 2 = p/X 
then eB ^ gb and the \Di%p\ 2 term in the Hamiltonian (|10|) leads to an infinite energy (the 
Meissner effect). If, on the other hand, we force eB = gb (which must then be spatially 
constant) then \tf)\ 2 = fp/X which is not at the minimum of the potential. 

The solution to the above problem is as follows: Suppose we first consider the config- 
uration described above in which \ip\ 2 = fp/X is not at the minimum of the potential. It 
may be surprising, but the fact is that this configuration is an extremum of the Hamiltonian 
(fLOl) subject to the constraint ([□]). To see this notice that this configuration does satisfy 
the equations of motion (||-0) but with 

gao = y/(l-/) (13) 

which is spatially constant. We see from Eq. (||) that the constant part of a$ behaves in the 
same manner as a shift in the chemical potential thus allowing the minimum of the potential 
to be consistent with the cancellation of eB and gb. 

It is also clear intuitively why this configuration is an extremum of the Hamiltonian 
(PH|). The main candidate for a fluctuation which might lower the energy is one in which 
ip is changed locally by an infinitesimal amount in the direction of the minimum of the 
potential. This will cause an infinitesimal amount of flux (J(eB — gb)) to thread through the 
system. If this flux is truly infinitesimal it will be less than one flux quantum and thus its 
gauge potential will lead to an infinite contribution to the gradient term in the Hamiltonian. 
(It is of course possible to consider fluctuations which have zero total flux but there is no 
reason to believe that these would lower the energy. The fact that they leave the energy 
unchanged to lowest order is evidenced by the previous argument that this configuration 
solves the equations of motion.) 

The configurations in which the flux is not infinitesimal but rather consists of an integral 
number of flux quanta is a finite energy excitation of this system. If, for example, we imagine 
modifying the chemical potential away from the special value we might hope that these 
vortices would lower the energy of the system and thus describe the Fractional Quantum 
Hall system away from special filling as a collection of these vortices. Unfortunately this 
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is not the case. We shall see shortly that for small deviations of /x from its "special" value 
these vortices always increase the energy. Let us first discuss why this is the case and then 
explain its consequences. 

The main point is that the configurations which are extrema of the Hamiltonian (|10D 



subject to the constraint ( |TTD are completely independent of the value of the chemical po- 
tential p\ We have, in fact, proven this already. Configurations are extrema of (|T0l) subject 
to ( JTT| ) if and only if they satisfy the equations of motion Thus if a particular ip is 

such an extremum (it may be a constant or, more generally, a multivortex configuration) 
with a given value pi of p then it will satisfy the equations of motion for some function . 
The same ip will satisfy the equations of motion for any other value P2 of p but this time 
with a new = clq + [1% — P\- It is thus also an extremum of the Hamiltonian with this 
new value /X2 of /x. 

Even though the extremal configurations for differing values of p are the same, the 
energetics may differ for different values of p. Notice, for example, that for p such that 
p/X = eB/(2nn) (let us call this value of the chemical potential p ), the configuration 
\ip\ 2 — A*oM clearly has the lowest possible energy. A single vortex configuration ip v (x) 
which solves the equations of motion (see |13|]) will have a larger energy than the ground 
state. Let e v be the excess energy of the vortex with respect to the ground state of the 
Hamiltonian with /x = /xo- Now consider an alternate Hamiltonian with /x = /xi ^ The 
same configuration ^ v (x) will still be an extremum of this Hamiltonian but its energy 
(which is the difference in energy between the vortex configuration and the configuration 
with \tp\ 2 = /x /A) will differ: 



e v ~ (a*i - A*o) J d 2 x \\i}) v {x)\ 2 - = e v - {p x - p )N v (14) 



where N v = ±l/n is the "particle number" of the vortex. Note that when /x is decreased it 
is preferable to form a vortex (N v < 0) whereas if p is increased an antivortex is preferred. 

Equation ( JH| ) has the following consequences. For small values of 5fi = /xq — A*i the 
vortex configuration increases the energy of the system at /x = We thus expect that the 
constant configuration with |-^| 2 = /Xo/A will be the configuration of lowest energy despite 
the fact that the potential energy is not at its minimum. This configuration has a density 
p = Po/X so that the system remains at the special filling fraction even after p has been 
shifted from p to p,\. (We emphasize again that this occurs for small shifts /x — p\.) It 
follows that in the Mean Field approximation 

T =° (15) 
dp lB 

for a range of p near p = p . This equation is familiar from the integer Quantum Hall 
effect and is due to the presence of a gap in the spectrum. (See Refs. [0] for a discussion 



of this with respect to the FQHE). It implies that a finite change in the chemical potential 
is required before the density can be modified. Returning to Eq. ( |T4j ) we see that when 
5p = \po — p\\ > e v /N v a single vortex has lower energy than the constant configuration 
\ip\ 2 = /xo/A. A weakly interacting gas of such vortices will have an even lower energy. 
We thus expect that near this value (e v /N v ) of 5p the lowest energy configuration of the 
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Hamiltonian (|10D subject to the constraint ( pTT]) will be a collection of vortices (which will 
likely form a lattice in the Mean Field approximation). 

The situation is similar when the magnetic field is modified instead of p. If we begin at 
the special filling fraction with eB = 27tnp /\ and change B at fixed p by a small amount, 
the lowest energy configuration of the Hamiltonian will occur at a new density for which 
2np/eB is still equal to 1/n but which will now not be at the minimum of the potential. 
As the magnetic field is increased (or decreased) further (again at fixed p = po) the cost 
in energy of a single vortex (or antivortex) becomes progressively smaller until at some 
critical value of the field it becomes negative. At that point the lowest energy mean field 
configuration is no longer a constant but rather a lattice of vortices in which case 2%p/eB is 
no longer equal to 1/n. (If the magnetic field is increased then the vortices will "condense" 
whereas if it is decreased the antivortices will "condense".) Thus if 2np/eB is plotted either 
as a function of /i or as a function of B there is a plateau surrounding the value po for which 
this ratio is constant and equal to 1/n. If, on the other hand, the B is varied at fixed density 
then we move off the plateau and the lowest energy mean field configuration consists of a 
vortex lattice. 

In the remainder of this paper we shall look more closely at the vortices of this model 
which, as we have discussed, will be solutions both for p = p and for values of p differing 
from p . We shall present some numerical solutions for these vortices which will allow us to 
estimate the value of p at which a lattice of vortices starts to form. The shape of the vortices 
will also lead to an estimate of the density at which the collection of vortices becomes non- 
dilute. The method we have chosen for finding vortices is by considering the Hamiltonian 
and constraint given by Eqs. (|I0|-|TT|) at a value of p = po = eB\/(27rn). We then look for 
radially symmetric configurations (which necessarily carry an integer number of flux quanta 
of eB — gb) which minimize the Hamiltonian. 

Anticipating the fact that our solution will be a vortex with an integer number of flux 
quanta we organize a radial ansatz as follows: First write 



ip(r, 6) = £(r)e~ 



ike 



(16) 



with k an integer. The Hamiltonian fllO| ) can now be written as: 



H 



2nrdr 




with the constraint 



— J fb(f)df\ £ 2 (r) 



-^ 2 W + ^ 4 (r) 



(17) 



r) 



2nn 



(18) 



(Note that if we were considering a value of p not equal to po we would still use the above 
equation but with p = po/f so that the Mean Field Solution would be £ 2 = fp/X.) We now 
define the function h(r) via the formula 



gb(r) — eB 



h'(r) 



(19) 
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where h'(r) = dh/dr. h(0) can be chosen equal to without loss of generality. The second 
term in Eq. ([L7|) is then proportional to: 




r'b{r')dr' 



h(r) 



(20) 



In order for the integral to be finite at large r we require h(oo) = k (or, more precisely, 
we require it to be an integer and we choose k in Eq. (|i~6|) to be that integer). (We only 
consider the cases k = ±1 in this paper, since those configurations have the lowest energy.) 
Furthermore, as is standard for all vortices, £ 2 (r) must vanish at the origin in order for the 
energy to be finite. The Chern-Simons condition Eq. ([18]) then implies 



27m 



27m 



as r — ► 



(21) 



This is the most difficult condition to implement in a numerical scheme in which the function 
h(r) is varied to minimize the energy. 

With the above definitions the Hamiltonian is given by 



H = 2tt / rdr 
Jo 




+ 



h(r) 



2m \ r 



e 2 (r)-^ 2 W + ^ 4 (r) 



(22) 



with h(r) chosen so that h(Q) = 0, h(oo) = k and £ 2 (r), defined by Eq. (PT|), is > 0. Notice 
that these conditions guarantee that the total flux of the vortex 



J d 2 x [gb - eB) = 2nk 



(23) 



The final step is to subtract, from the energy of the vortex solution of Eq. (|22|) , the energy 
of the Mean Field solution £ 2 = ffi/X. This results in a vortex energy given by: 



H v = 2tt 



POO 




/ rdr 


±1 


Jo 


2m \ 




"Ml - /) ^(r) -f£)+± (e(r) - f 




(24) 



The procedure at this stage is to search, numerically, through the space of such functions 
h(r) until the Hamiltonian is minimized. 

One point which is clear is that the form of the vortex solution (for which J (eB — gb) > 
is quite different from that of the antivortex solution. The reason is that £(r) and thus b(r) 
must vanish at the origin. As a consequence the density ^ 2 (r) cx b(r) for the vortex solution 
can be a monotonic function of r which increases from at the origin and reached eB at 
infinity. The antivortex solution must however be zero at the origin then increase to a value 
greater than eB (so that J(eB — gb) < 0) and then decrease again to attain its asymptotic 
value eB as r — > oo. We shall of course see this behavior clearly in the numerical solutions 
shown in Figures 1-4 below. 

For the numerical work we chose some representative values for the parameters of the 
model M: 
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p = .010 eV 
Po/X = 10 3 pm~ 
m = .08 m e 



-2 



(25) 



(In the units K — c — 1 this translates to A = .00025 eV" 1 , p = /io/A = 40 eV 2 , eB = 
250neV 2 and m = 41000 eV.) Note that po/X is the density of carriers. (Although the 
above values were chosen to be representative of experiments which exhibit the FQHE it 
is difficult to call them "realistic" since the CSLG model has been greatly simplified by 
replacing the Coulomb interaction of the charge carriers with a short range interaction.) 

In Fig. [TJ and Fig. |2] we present the numerical solution for the vortex configurations (in 
which gb is lowered relative to eB and h(r) — > — 1 as r — > oo). We plot the functions h(r) 
and the density p{r) for v = 1/3, 1/5 and 1/7 respectively. When changing z/, the density p 
remains fixed as the magnetic field B is varied. Figures |3| and |] contain plots of h(r) and 
p{r) for the same values of v but now for the antivortex configurations. In Table | we 
present the energies and some measure ro of the size of each vortex and antivortex. We have 
arbitrarily chosen the size of the vortex as the value of r at which the energy density has 
reached 99% of its total value. 

We are now ready to describe quantitatively (within this model) what happens when the 
chemical potential is varied from pq. As discussed in great detail in this paper there is no 
change in the density unless p — po is approximately equal to the energy of a vortex times n 
(i.e. the energy per particle of the vortex). We can now see, quantitatively, how this works 
from Eq. ( p4| ) . If / ^ 1 so that p = po / f ^ po then the energy of the vortex is simply 



where the minus sign is for a vortex and the plus sign for an antivortex. Thus for p < p the 
vortex configuration has lower energy than the antivortex configuration. We naively expect 
that when po — p = ne v (po) (or near this point) the Mean Field ground state should be 
a condensate (possibly a lattice) of vortices. Conversely when p > p the antivortex has 
lower energy and when p — p = ne v (p ) we naively expect a condensate of antivortices. 
Unfortunately, for our choice of parameters, e v is quite large. Thus the value p cr or p at 
which this condensate occurs in the above naive calculation and which is shown in Table | 
differs from po by an unreasonably large amount. This leads, in particular, to a negative 
value of p c r for the n = 5 and n = 7 vortex configurations. In fact, as p is varied from po 
towards p cr , another critical value p of p is reached at which 2-apjXeB = l/(n + 2) well 
before p cr is reached. At p the system is better described by a Chern-Simons theory with 
the new value v — l/(n + 2) of the filling fraction. 

In light of the above remarks we should try to understand whether in fact one does form 
a vortex condensate in our model at our chosen values of the parameters. Certainly the 
vortices must be strongly interacting (or overlapping) well before p or p cr is reached. This 
can be better understood by first supposing that such a condensate is formed as p is lowered. 
The approximate density p\ at which a description of this condensate in terms of the single 
vortex solutions presented above fails depends most prominently on the size of a vortex. 
In Table we show the approximate density of vortices p v \ and chemical potential p\ at 
which the vortices begin to "touch" . For a hexagonal lattice of vortices this will occur when 




(26) 
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p vl = l/(2y / 3<i 2 ), ie, where d (the distance between vortices) equals the size of the vortices. 
Table also shows the corresponding values of the filling fraction. Notice that in most cases 
the vortex condensate becomes dense well the before the "next" value of n (i.e. well before 
I'npjeB = l/(n ± 2)). We conclude from this that an approximation in terms of a dilute 
gas of vortices breaks down well before /i = fi cr . It is thus likely that even for our chosen 
values of the parameters a condensate of vortices will form in the Mean Field description. 
The formation of this vortex condensate and the resultant pinning of the vortices is what 
gives rise to the hall plateaus by allowing the system to continue to behave as if it were in 
a state v — 1/n even after the magnetic field or the chemical potential has been changed to 
move it away from that value. This is because the (anti)vortices accommodate the localized 
excess (deficit) of charge. 



Summary 

In this paper we have studied the Mean Field behavior of the CSLG description of 
the FQHE when the filling fraction deviates form the "special" filling fraction for which 
v — 1/n with n and odd integer. We have shown how the Field Theoretic description of 
this model at a fixed chemical potential /i and magnetic field B can be studied for a range 
of /j, surrounding the value fiQ corresponding to the special filling fraction. For small values 
of \fi — {Iq\ (and at zero temperature) the density is independent of fi. This is reminiscent 
of what occurs for the integer Quantum Hall Effect. As /i is decreased beyond some /i~ 
we show how the homogeneous Mean Field configuration is unstable to the formation of a 
condensate of vortices. If /x is increased above some /z+ the instability is to the formation 
of a condensate of antivortices. We have presented a numerical example of these vortex and 
antivortex configurations and we estimated the densities and filling fractions at which the 
description in terms of a noninteracting system of vortices breaks down. 
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TABLES 



vortex/antivortex 


n 


size (fj,m) 


energy (xlO 2 eV) 


Ha- (xl0~ 2 ey) 




3 


0.027 


0.25 


0.25 


vortex 


5 


0.018 


0.21 


-0.05 




7 


0.015 


0.20 


-0.4 




3 


0.035 


0.83 


3.5 


antivortex 


5 


0.025 


0.80 


5.0 




7 


0.020 


0.79 


6.5 



TABLE I. Energy and size (ro) of vortices and antivortices for n = 3, 5 and 7. \i cr is a naive 
estimate of the value of the chemical potential at which a condensate of these configurations is 
expected to form (/xq = W~ 2 ev). fj, c r is more carefully described in the text. 



vortex/antivortex 


n 








3 


400 


2.6 


vortex 


5 


890 


4.1 




7 


1280 


5.7 




3 


240 


3.2 


antivortex 


5 


460 


5.5 




7 


720 


7.7 



TABLE II. The density p v \ and the corresponding filling fraction v\ at which the condensate 
of vortices is expected to become dense. (See text for a precise definition.) 
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FIGURES 
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FIG. 1. The function h(r) corresponding to a vortex for u = 1/3, 1/5 and 1/7. 
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FIG. 2. Density profile of the vortex configuration for v = 1/3, 1/5 and 1/7. 
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FIG. 3. The function h(r) corresponding to an antivortex for v = 1/3, 1/5 and 1/7. 
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FIG. 4. Density profile of the antivortex configuration for v = 1/3, 1/5 and 1/7. 
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